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We present the first 3+1-dimensional numerical simulations of scalar fields with non-minimal 
kinetic terms. As an example, we examine the existence and stability of preheating in the presence 
of a Dirac-Born-Infeld (DBI) inflaton coupled to a canonical matter field. The simulations represent 
the full non-linear theory in the presence of an expanding Universe. We show that parametric 
resonance in the matter field, along with self-resonance in the inflaton, repopulate the Universe 
with matter particles as efficiently as in traditional preheating. 



Scalar fields need not be Klein-Gordon in four dimen- 
sions; theories that include extra dimensions, whether 
large or small, have effective four-dimensional kinetic 
terms that are non-minimal. This has led to an increased 
interest in the dynamics of non-minimal models both to 
describe inflation, e.g. Dirac-Born-Infeld (DBI) infla- 
tion [ij and, even more recently, dark energy In 
this paper we will focus on DBI inflation, as an example 
of a system for which we have a scalar degree of freedom 
whose behavior is self-consistent and stable, and whose 
non-minimal behavior is central to the dynamics of the 
model. Here we focus on what happens at the end of in- 
flation when the scalar inflaton becomes inhomogeneous 
and its couplings to other fields are important. To our 
knowledge, these are the first three-dimensional lattice 
simulations of this type of field theory. 

Preheating [5-17| provides a mechanism by which the 
cold post-inflationary Universe can quickly and efficiently 
transfer energy into a matter sector, via a period of para- 
metric resonance or a regime of tachyonic instability. For 
the most part, studies of preheating have focused on the 
existence and stability of these processes in the presence 
of different inflationary potentials, e.g. llTI- 22] . multiple 
flelds [23], and multiple decay channels 12J]. DBI infla- 
tion might end with a coherently oscillating scalar fleld, 
and it has been unclear whether preheating can persist in 
the presence of non-linear terms in its equation of motion. 
The search for preheating in non-minimal models began 
in [i^ , where the authors showed that a canonical scalar 
fleld can enhance the effects of parametric resonance in 
a coupled non-minimal matter fleld. 

The flrst studies that considered non-minimal infla- 
tionary flelds [2^ looked at the production of particles 
during the flrst oscillation, instant preheating, as a pos- 
sible source of non-gaussianities in the Cosmic Microwave 
Background (CMB); restricting the treatment to the first 
oscillations allows a complete analytic treatment of the 
amplification of the matter field, but does not cap ture the 
majority of the preheating dynamics. Later, j27| studied 
a DBI-type model with perturbative departures from a 
canonical kinetic term, so that the potential dominated 
the preheating dynamics. In this case preheating per- 



sists and is perturbatively enhanced, but in addition self- 
resonance occurs in the inflaton. At the other extreme, 
28| studied a case where the non-minimal kinetic term 
is dominant. In this case the potential is irrelevant dur- 
ing most of the oscillation and the fleld proflle develops 
a "saw-tooth" oscillatory proflle. Further, the authors 
argue that this feature suppresses preheating. 

Here we study the full non-linear theory on a three- 
dimensional lattice. We will not assume that the non- 
minimal terms are small, nor that inflaton is purely ho- 
mogenous. Such a treatment allows us to study the exis- 
tence of resonance in the coupled fleld, any self-resonance 
in the inflaton, and any non-linear effects that take place 
after the resonant periods cease. Further, all analytical 
analyses of preheating with non-minimal kinetic terms 
require many simplifying assumptions: they may disre- 
gard Hubble friction and must ignore the spatial gradient 
terms in the equation of motion. Numerical simulations 
are needed to verify the possibility of preheating in this 
model outside the non-relativistic limit, and to see any 
interesting physics that arises from the interactions be- 
tween the two flelds. We begin by considering the stan- 
dard action for a DBI inflaton [H 
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We will take the warp factor / to satisfy [l[ 
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It is related to the geometry of the extradimensional 
space in which the brane is moving; the form in Eq. [T] 
corresponds to the case of a cut-off AdS throat. In a 
braneworld scenario would correspond to the radial dis- 
tance of the brane to the horizon. Note that, in the limit 
where f X <^ 1, we recover the canonical kinetic term; the 
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magnitude \fX\ therefore characterizes departures from 
the standard scenario. 

The full non-linear equation of motion for 0, from 
Eq. m is 
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in a Friedmann-Lemaitre-Robertson- Walker, 
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background. The combination 7 = (1 — fX) ' de- 
scribes how relativistically the brane is moving and is 
analogous to the usual Lorentz factor. We will assume 
that the matter field is a massless canonical, Klein- 
Gordon field whose equation of motion is 
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and the full potential of the model is the combination 
of the inflationary potential and a coupling between the 
two fields, 



T/(</),x) = imV + ^.gVx'- 



(7) 



Since we arc interested in effects due to the non-minimal 
nature of the field, we will not allow the parameters in the 
potential to change from the standard values used in the 
preheating literature, m = 10~^ uipi and = 2.5 x 10~^. 
We will take the ratio X/n'^ as a free parameter. 

The contribution to the energy density from the infla- 
ton. 
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differs from the energy density associated with the matter 
field. 
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The combination of the two with the potential energy, 
p = P4, + Px + V{(j), x)i together source the self-consistent 
expansion of the Universe, 
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(10) 



Numerical Simulations: Simulations were performed 
using Grid and Bubble Evolver (GABE) [111, a new lat- 
tice evolution program that uses a second-order Runge- 
Kutta metho d of integration. Unlike existing lattice 
codes 



33j that use symplectic integration routines, 
GABE stores the field and field derivative value at the 
same times during the time step. This method is slower 
and requires more physical memory to run than do sym- 
plectic integrators but is necessary when terms include 
the product of the field with its time derivative. We offset 
the slower performance by parallelization using OpenMP. 
We postpone a complete analysis of all relevant values of 
X/ p'^ for a subsequent publication, instead focusing on a 
single case, p = lO'^mpi and X/p^ = 5000 mpj**, as a rep- 
resentative model. In the numerical simulations here, we 
use a box of = 256'^ points. Running on four cores, 
each run takes approximately 8 days to complete. 

We initialize the homogeneous mode of the matter field 
as x{k = 0, i = 0) = 0, whereas the homogeneous modes 
of the inflaton and its derivative are consistent with their 
values at the end of inflation, when a drops below zero. 
We calculate these initial conditions using a sixth-order 
Runge-Kutta homogeneous field evolution program to 
evolve the field, (j>, during infiation, when it obeys Eq. |3] 
without the gradient terms. For our case of interest here. 
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the dynamics of the field is different, even at the end of 
inflation, these values are close to those expected in the 
canonical model. 

Both fields and their derivatives are initialized with 
fluctuations that are consistent with the Bunch-Davies 
vacuum, in a way analogous to (soj . so that the power 
spectrum of the fields on the initial slice are 
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After initialization, the fields are evolved for several hun- 
dred oscillations of the inflaton fleld, until tf ^ 250 m~^, 
in an expanding background. 

The software calculates the means, variances, and en- 
ergy density of both fields as well as their energy and 
power spectra. Two-dimensional slices of the field profile 
arc also generated, as well as the space-averaged value of 
7. The scale factor grows by a factor of about Uf ^ 17 
(when X/ p^ ~ 5000) compared to ay ~ 27 (in the canon- 
ical case for the same final time). 

Results: The evolution of the inflaton field during both 
infiation and its first few oscillations show the effects ex- 
pected due to the presence of these non-canonical terms. 
In the DBI the velocity of the field grows, the rel- 

ativistic factor 7 also grows, decreasing the infiuence of 
the potential term on the acceleration of the field. Math- 
ematically, to guarantee the reality of 7, cannot exceed 
j-1/2. iiowever, in practice we see that the true speed 
limit is somewhat smaller and is a consequence of the 
non-linear dynamics of the field. We can easily see the 
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effects of the speed limit on the oscillations of the infla- 
ton after inflation. In fact, the mean value of (j) on the 
lattice does not oscillate sinusoidally, but instead reaches 
its maximum velocity quickly, after which it travels at a 
constant speed for much of each oscillation. The evolu- 
tion of (p thus approaches a saw-tooth pattern as A//x^ 
increases, as predicted in [2I]; see Fig. [T| 
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FIG. 1: The evolution of the mean inflaton field value for 
simulations with X/j-i'^ = 5000, 70 ~ 8.7 (red, solid), and the 
canonical case, / ^ 0, (blue, dotted). 

We have run simulations with initial values of 7 varying 
from 1 < 7 < 10, consistent with the Planck 2013 con- 
straint 7 < 14 at 95% confidence [s^l- In all simulations, 
we see significant effects due to parametric resonance by 
approximately t 100 m~^. We can identify parametric 
resonance by the exponential amplification of particular 
modes of the matter field, x, which result in an exponen- 
tial increase in the variance of the matter field over time. 
We can see this schematically by noting the inflaton is a 
coherently osciUating field, cj) = ^{t)', the mode equations 
for the matter field, 

XS + ^lxj:+{k'+g'^')Xi:^0, (12) 

are then damped harmonic oscillators with a time- 
dependent mass. In the case of a sinusoidally varying $ 
we can reduce Eq. 1121 after ignoring the expansion of the 
Universe, to the Mathieu equation and predict the spec- 
trum of amplifications. If we allow $ to be a "saw-tooth" 
function whose amplitude decreases and whose period in- 
creases, the consequences for preheating are unclear. On 
one hand, we expect that the time-varying period of os- 
cillation should do harm to the period of parametric reso- 
nance. Some modes will experience small amplifications 
during each oscillation, but there is no assurance that 
any particular mode is amplified repeatedly. 

On the other hand, the "saw-tooth" is actually many 
Fourier modes; it represents many forcing terms, each 
with different resonance bands. Since the resonance in 
this case is much broader, we might expect that more 



modes are in resonance at any given time and the effi- 
ciency of parametric resonance is increased. Fig. [2] shows 
the comparison between the canonical and the DBI cases. 
The X field is efficiently amplified during the preheating 
process, not hindered by the "saw-tooth" oscillations. 
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FIG. 2: The variance of the matter, x, fields. The red (solid) 
lines represent the fields when A//i^ = 5000 (initial 7 = 8.7), 
and the blue (dotted) lines represent the canonical (/ — >■ 0) 
case. 

Perhaps more interesting is the burgeoning importance 
of self-resonance. The extra terms in Eq. 2] give rise 
to self-interactions that provide a mechanism for self- 
resonance. Unlike in the canonical case, the modes of 
(j) undergo strong self-resonance in the presence of non- 
minimal terms. Fig. [3] shows the dramatic difference in 
the two regimes; indeed, for the case of interest here, 
we see that the self-resonance is faster and more effi- 
cient than the induced parametric resonance in the mat- 
ter field. This can be seen in Fig. |31 since the variance 
of (j) for the non-minimal case grows quickly in the early 
stages of the simulation, whereas in the canonical case, 
the variance of (j) decreases during this time. The exis- 
tence of self-resonance is generic for different values of 
A///'*, becoming more significant as we depart further 
from the canonical case. 

In both simulations, the periods of resonance cease 
around t ^ 25 m~^. At this point, the two fields are suf- 
ficiently inhomogcneous that they begin to interact and 
thermalize via the g^4P"x^ /2 interaction. Lack of three- 
leg interactions, as is common in preheating simulations, 
makes it impossible for the (f> field to decay completely 
Nonetheless, we see a final state consistent with any 
simulations of preheating after quadratic inflation. 

Discussion: The above simulations, the first to model 
a scalar field obeying an equation of motion with non- 
minimal kinetic terms in a 3-|-l-dimensional Universe, re- 
veal that even when those non-canonical terms are large, 
preheating occurs. We see that the parametric amplifi- 
cation of non-thermal modes of a coupled matter field 
are just as efficient as in the traditional preheating sce- 
nario. Moreover, as we depart from the canonical case. 
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FIG. 3: The variance of the inflaton, 0, field. The red (sofid) 
fines represent the fields when \/ = 5000 (initial 7 = 8.7), 
while the blue (dotted) ones represent the canonical (/ — 0) 



parametric resonance is no longer the primary cause of 
preheating; instead, the self-interaction of the field causes 
almost immediate and efficient particle production. The 
fact that the inflaton does not vary sinusoidally after in- 
flation is not a death sentence for preheating; on the 
contrary, the spectral diversity of the homogeneously os- 
cillating mode of (j) acts as a source of inhomogeneity 
equally as efficient as a pure sinusoid. 

Since the structure of preheating is so similar, one can 
additionally expect that the gravitational wave spectrum 



from preheating |35l - l43| . as well as any other observable 
consequences, e.g. non-gaussianities 44 "i^, should not 
be any less significant, although we leave these analy- 
ses for other work. We also delay a discussion of how 
self-resonance effects are important when considering the 
generation of pseudo-stable, non-linear structures in the 
inflaton field ^HH. 

Just as importantly, this work serves as a proof-of- 
concept that 3-|-l-dimensional lattice simulations can 
simulate fields with non-minimal kinetic terms. This 
should pave the way for simulating tensor fields with non- 
minimal kinetic terms sourced by scalar fields. 

In a future paper, we will present the results of simu- 
lations with a range of values of 7 at the end of inflation. 
We also delay predictions for the spectrum of amplifica- 
tions that one expects from this process. 
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